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SUMMARY 


A method is given for the design of three-dimensional nozzles uti- 
lizing axisymmetric flow. The nozzle can be designed to produce imlfoiin 
flow in a test chamber of arbitrary cross section. The method, is applied 
to obtain the final coordinates of a Mach number 10 nozzle for which a 
square test section is specified to reduce the possibility of axisym- 
metric Imperfections at the wall and to provide for the installation of 
schlieren windows . Radial flow is used in a portion of the flow field 
to reduce the computation time. The remainder of the flow field is 
computed by the method of characteristics, but a simplified method is 
used near the axis. Tables which facilitate computation of the radial 
flow and the flow near the axis are included. Transition streamlines 
determined from the analytic expressions of Kuno Foelsch are compare’d 
with the streamlines obtained from the characteristics net of the Mach 
number 10 nozzle. The Foelsch streamlines deviate from the flow-net 
streamlines by as much as 12 percent. Similar analytic expressions are 
derived from the geometric properties* of the flow. These new expressions 
result in transition streamlines with a maximum error of about 4 percent. 


INTRODUCTION 


Recent research has indicated that three-dimensional supersonic 
nozzles may become more desirable for high Jfech number tunnels than 
conventional two-dimensional nozzles. For example, in two-dimensional 
nozzles designed for test Mach numbers much greater than 5^ the flow is 
very sensitive to any change of the extremely small dimensions at the 
minimum section. The high temperatures required to avoid liquefaction 
of the air at these high Mach numbers make the problem of obtaining 
dimensional stability of the small slit-like minimum extremely difficult. 
(See reference 1.) In addition, the excess growth of boundary layer 
along the center of the nozzle side plates may also interfere with the 
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design flow. A two-stage nozzle discussed in reference 2 avoids the 
first difficulty hut fails to operate satisfactorily primarily because 
of the excess boundary layer. Consideration of three-dimensional nozzles 
thus becomes imperative not only for tests of stationary models at the 
higher Mach numbers but also for ballistic tests where the projectile 
must pass through the minimum of a supersonic nozzle. 

A method for the design of three-dimensional nozzles based, on axi- 
symmetric flow is presented in this paper. Although the flow is axi- 
symmetric throughout, the final cross-sectional shape of the nozzle may 
be arbitrary. The general concept for the determination of the arbi- 
trary cross-sectional shape was applied to the design of an axisymmetric 
nozzle at the Langley Aeronautical Laboratory by Mr. Morton Cooper in 
19^7- The design method presented in this paper is general; however, 
as an illustrative example the design of a Mach number 10 nozzle with a 
square test section is included. Specification of a square cross- 
sectional shape for the test chamber provides for the installation of 
conventional schlieren windows, A further advantage of the resulting 
noncircular walls is that the possibility of incurring axisymmetric 
in^jer feet ions at the nozzle wall is reduced. Disturbances originating 
from axisymmetric imperfections focvis along the axis of the nozzle. 

This focusing effect has been known to cause disturbances of considerable 
magnitude in the vicinity of the nozzle axis. A theoretical analysis of 
the three-dimensional focusing effect is given in reference 3^ part II. 


LIST OF SYMBOLS 


A,A',B 
C,D,E 7 
P,K 

c 

M 

V 


x,R,0 

s 


points in axisymmetric flow field as indicated in figure 1 
or sketch on p. 10; used as subscripts for conditions at 
these points 

local speed of sound 

Mach n\imber 

velocity (with bar to indicate vector quantity) 

distance measured from source point in radial flow field 

distance measured from the source point to sonic sphere in 
radial flow 

cylindrical coordinates for flow field and nozzle 

distance along a Mach line measured in meridian plane from 
the axis of symmetry 
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x,y,z 

7 

e 

P 

i 


Cartesian coordinates for nozzle 
Mach angle, sin“^ I/m. 

ratio of specific heats = 1.400 throughout 
local flow angle 

total expansion angle integrated from the sonic line 
in meridian plane 

density- 

stream function for a compressihle axisymmetric flow 
dimensionless stream function 


Subscripts ; 


r radial flow 

o isentropic deceleration to zero velocity 


NOZZIE DESIGN 


For convenience in discussion, the nozzle is divided into four 
regions: I, approach; throat; III, radial flow; and IV, transition 

region, as indicated in figure 1. 

Approach . - The specification of parallel flow at the mini m um 
section facilitates the computation of the supersonic portion of the 
flow. Theoretical methods (references 4 and 5) are available for the 
design of an approach to give parallel flow; however, it is usually 
more conservative to base the approach- sect ion design on existing wind 
tvmnels which have the desired characteristics. The coordinates for 
the approach of the Mach number 10 nozzle were therefore adapted from 
those used in the Langley 8-foot hi^-speed tunnel. In the absence of 
a better criterion the square approach section was obtained by using 
the same rate of area change -with axial distance as in the circular, 
approach of the Langley 8-foot high-speed tunnel. 

Kadlal flow.- Construction of the supersonic flow field may start 
with either an arbitrary streamline A'C (fig. l) together with an 
assumed flow distribution in the minimum or with the specification of 
a flow distribution along the axis of symmetry of the flow. The second 
specification was chosen in this investigation primarily to expedite the 
computation procedure; however, if this center— line distribution were 
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entirely arbitrary, the extent of a physically possible flow would he 
severely limited; that is, Mach lines of the same family would converge 
too rapidly or intersect before the computation had proceeded far 
enough fj’om the axis to result in a we3-l -proportioned nozzle. This 
difficulty is avoided by specifying radial flow along a portion of the 
axis such as line BD. If this specification is used, the flow can be 
immediately determined for region III from the radial-flow equations 
and, in general, a physically possible flow results for region II. 

The use of radial flow in a part of the flow field has the addi- 
tional advantage of reducing the computing time, since the flow proper- 
ties along the bounding Mach lines BC and CD can be determined exactly 
without recourse to the step-by-step method of characteristics. The 
equations for three-dimensional radial flow .are based on the assumption 
of spherical symmetry. The resulting simplification of the general 
equations of motion allows their integration in closed form as indicated 
in reference 6 and 7- In addition, it is possible to integrate the 
characteristic equations ae in two-dimensional steady isentropic flow. 
For convenience, the three-dimensional radial-flow equations based on 
the general steady-flow relations (for example, reference 8) are pre- 
sented in the Appendix. Values have been calculated by use of these 
radial-flow equations and the results are tabulated for a range of 
parameters from M = 1.0 to M = 11.93 in table I. This change in 
M corresponds to a change in expansion angle 0j, which is used as the 

argument, from 0° to 53.375°. 

If point C is at the intersection of the Mach lines bounding the 
radial- flow region, as illustrated in figure 1, it will be the inflec- 
tion point on the streamline A'CE. This streamline may be taken as 
the wall of a nozzle, and the maximum weQ.1 angle and the coordinates of 
the Inflection point can then be determined immediately from the local 
flow angle at C. The value of the flow angle at C is entirely arbi- 
trary; however, a reasonable value can be determined as a compromise 
between the requirements for over-all nozzle length and maximum rate of 
expansion. For the Mach number 10 nozzle, the flow angle at C was 
taken as l6°. 

Throat .- The flow in region II (the, throat) is computed by the 
method of characteristics with the initial conditions of a known flow 
distribution along line BC from radial flow and .an arbitrary distribu- 
tion along line AB as shown in figure 1. The Mach number 10 nozzle 
computation was simplified by taking the arbitrary distribution as a 
straight line with a slope given by the radial flow at point B. In 
general, this assumption is incompatible with a straight sonic line; 

/ -siA 

indeed, a necessary condition for a straight sonic line is ( r— ) = 0. 

\ox/M=l 

However, for high Mach number nozzles, this simplifying assumption 
results in a satisfactory solution because of the small dimensions of the 
minlTtmm section. The characteristics net for region II is shown in fig- 
ure 2 with a streamline throtigh point C included. 
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The characteristic equations may be used in a form where the more 
usual velocity parameter, such as that used in reference 9, is replaced 
by the Mach number by using the differential relation between Mach 
number and velocity. Sample calculations indicate that the use of this 
parameter for the high Mach number range reduces the numerical error 
inherent in the method of characteristics. The particular parameter 
used makes little difference at the lower Mach numbers. 


T'/hen the Mach number is used, as a parameter/ the characteristic 
equations are 


da + — dM t -j— : . ■ } dx = 0 (1) 

+ - 1 cot a 7 1 ) 

and 

dR = tan(a t p) dx (2) 


where the upper signs correspond, to the left Mach' lines and the lower 
signs to the ri^t Mach lines. The distinction between these Mach lines 
is illustrated in figure 2. The last term on the left-hand side of 
equation (l) becomes indeterminate at the axis of symmetry, that is, 
as R and 6 approach zero. If the limit of this term is taken 
(reference 9 ) as R-»0, equation (l) becomes 


da + 




- 1 


Ml + 


7-1 




dM t 


1 - 1 
2 M^l + 2L. - . 1 1^2^ 


M dx = 0 

,^^/R=0 


( 3 ) 


Thus, if I — ] is specified along the center line, equation ( 3 ) may 
\^/R=0 

be used in conjunction -with equations ( l) and (2) to compute the flow. 

An alternative method of computing the flow in the vicinity of the 
axis can be used. This method consists in reducing equation ( 3 ) by 
dividing this equation by dx and using a relation which follows from 
the general expression for the total derivative of M with respect to x 
along a Mach line; that is, dM/dx = Bm/5x + (Sm/Sr) ( dR/dx) . If this 
total deri-vative is evaluated at R = 0, there results the required rela- 
tion, dM/dx = since in axisymmetric isentropic flow 

(Bm/Br)j^^ = 0. (See reference 3, part I.) 
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Hence the characteristic equations evaluated at the axis hecome 


d0 


+ 1 ^ dM = 0 

^ M^l + ^ ~ 


(4) 


Then equations (2) and (4) are used to compute the flow in the immediate 
vicinity of the axis; that is, at points 2,1; 3>2; 4,3; . . . , in fig- 
ure 2. Equation (4) is the differential form of the corresponding equa- 
tion given hy Goldstein as a finite-difference equation (reference 3, 
part I, page 2l8) . 


A sinrplificatibn is possible since equation (4) is of the same 
form as the corresponding equation for radial flow. (See Appendix.) 
Hence, the integral of equation (4) or the corresponding tabulated 
results (table l) can be used, in effect, to determine the hodograph in 
the region of the axis for any axisymmetric, isentropic flow. The 
method should be sufficiently accurate even at a small distance from 
the axis if 8 m/ 8R is smaJ.1. 


The row of points adjacent to the axis in region II (figure 2) have 
been computed by both equations (3) and (4) for the Mach number. 10 nozzle. 
The results of this computation are compared in table II. The differences 
between the two methods are small for all the points and tend to decrease 
as point 2,1 is approached. The remainder of the flow field in region II 
was then computed with equations (l) and (2) for both methods at the 
axis. A relative compaxison of the effect of the two methods on the rest 
of the flow field may then be obtained by comparing corresponding stream- 
lines through the two flow fields. If these streamlines are determined 
by the flow-angle method (described in a subsequent section entitled 
"Streamlines") starting from point C, for example, and 'working toward A', 
accumulative errors arise which prevent the streamlines from passing 
through the theoretically exact point A' . In this connection it is noted 
that the R coordinate of point A' can be determined exactly from the 
mass-flow requirement at, the minimum (assuming parallel flow at M = 1 
along a line through A’ normal to the axis); whereas the x coordinate 
of A’ depends on the assumption made for the flow distribution along the 
line AB. Thus, if the streamlines in both cases are started at the fixed 
initial point C, the error in s'b ^ 0.96 percent for 

the first method (equations (l), (2), and (3)); whereas for the alternate 
method (equations (l), (2), and the integral of equation .(4) ) , the error 
was 1.02 percent of the known radius at x = x^. 
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Transition . - The flow in region IV (transition region) can he com- 
puted by equations (l) and (2) from the initial conditions of a radial 
distribution along the Mach line CD and parallel and uniform flow along 
the Mach line DE. The characteristics net for region IV in the Mach 
number 10 nozzle is shown in figure 3 with a streamline through point C 
included. 

Mesh size.- The axisymmetric flow fields in regions II and IV 
(flg.TT were computed for the Mach number 10 nozzle by using the method 
of finite differences applied to equations (l) to C3)» The flow for 
region II was computed with seven steps along line AB and eight steps 
along line BC as shown in figure 2; the nearest point to the sonic line 
had a Mach number of 1.019928. Region IV was computed with 2° steps 
in 0j along CD and 10 uniformiy spaced steps along DE as shown in fig- 
ure 3- Three iterations were made for each point according to the pro- 
cedure of reference 9 so that the parameters were allowed to converge 
to nearly constant values. Altogether, this computation totaled 
7^ points and was performed on an electrically operated desk-type 
computing machine at an average rate of 3 hours per point for a total 
of approximately 220 hoxirs of computing time. Application of the 
simplified method (integral of equation (4)) at the axis would reduce 
the total conqjuting time to about 200 hours. Without the specification 
of radial flow along BD, 21 more points would have been required, and 
the total computing time would have been increased to about 290 hours. 

When a streamline through this flow net was determined, appreciable 
errors restilted because of the size of the steps chosen for the compu- 
tation. For this reason, computations were made with a finer flow net 
with steps approximately I/3 the size of those in the coarse net 
previously described and requiring a total of 573 points. The Bell 
Telephone Laboratories X-66f44 relay computer was used to perform this 
task. The available Bell computing tapes utilized the characteristic 
equations in the form given in reference 10 and made only two iterations 
for each point. The Bell computer performed the computation at a rate 
of 0.3 hour per point for a total of 172 hours. 

Sample computations have demonstrated that the values of the flow 
parameters will converge in the iteration process more rapidly if equa- 
tions (1) and (2) and the Iteration procedure of reference 9 ai’s used 
instead of the equations and iteration procedure of reference 10. How- 
ever, a comparison has been made of the accuracy of the flow nets as 
determined on the desk conq)uter and the Bell computer. The results of 
this comparison are given in table III, in wHich the ordinates of the 
streamlines determined by the flow-angle method are compared with the 
theoretically exact values at points C and A'. The errors result from 
the approximations inherent in the finite-step method employed in 
computing the flow net and are some fimctlon of the size of the steps. 
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Since the results obtained with the Bell con^niter were more 
accurate j in spite of the fewer number of iterations made and the slower 
convergence of the iteration procedure, the advantage of maintaining a 
fine net is obvious. For these calculations, the increase in accuracy 
is of the same order as the decrease in mesh size. 

Streamlines . - After the flow field has been computed the stream- 
lines may be obtained either from the local flow angles or from the 
stream function. 

Stokes' stream function I'’ (reference 8) for an axisymmetric 
compressible flow can be defined by the differential relation 



(5) 


where ds is the differential distance along a Mach line in the meridian 

plane. Introducing a dimensionless stream function i ^ into 

""cr Po<=o . 

equation (5) and integrating from the axis / = 0^ to gives the 


value of the stream function at 


•cr 


•cr 


■cr 


fs = 


'V^cr R 




( 6 ) 


where the isentropic flow relations for p/pq and c/cq have been 
used. Within the radial flow region BCD (fig. l) , t can be defined 
by the equation 


d.t o 

— = pVr^ sin 9 

d9 


( 7 ) 


Integrating along the arc and dividing by 

7-H 




( 8 ) 



KA.CA TN 2711 


9 


since p, V, and r are constants along the arc of radius r. Now 
the value of $ can be found for a given value of 9 along the Mach 
lines CD or BC from equation (8). Then, since $ is constant along a 
streamline, graphical integration of equation (6) along the Mach lines 
in region II or IV up to the desired value of f locates the remaining 
points for this streamline. In practice, the most accurate procedure 
is to obtain a value of $ for each point in the Mach net (that is, 
intersection points of ri^t and left Mach lines ) by integrating equa- 
tion (6) along the left Mach lines. Any desired value of $ and thb ^ 
corresponding x and R coordinates are then easily obtained by linear 
interpolation between the points. Integration along the left Mach lines 
is more accurate in region IV than integration along the ri^t lines 
because of the more gradual change in flow variables along the left lines 

The streamlines may also be obtained from the flow angles by 
extending the streamlines directly through the R - x plane by a step- 
by-step iteration procedure. This iteration procedure consists of 
extending the tangent to a streamline from a known point in the flow, 
such as a point on the Mach lines BC, CD, or DE, until another Mach 
line is intersected. As an exaniple, consider the streamline through 
point E, in figure 3- By linear interpolation between the known values 
of the flow angle at points 1 and a value (say Sp) obtained for 
this first intersection of the tangent line and the Mach line. The mean 
of 02 initial value 0g is then used as the slope for a new 

line through the initial point'. The intersection of this line with the 
Mach line 1-3 gives a new flow angle 02'» This process is repeated 
until the change in ©2 for successive iterations is no longer signi- 
ficant, whereupon the entire procedure is repeated in the adjoining mesh. 

The flow-angle method was not used for the final computation of the 
streamlines of the Mach number 10 nozzle because of the accxomulatlve 
errors at the end points A' and C. Nevertheless, the flow-angle method 
was found useful for comparing the effect of mesh size on the accuracy of 
the flow field. The stream-function method results in streamlines with 
the exact coordinates at points A', O, aud D. Furthermore, this method 
requires less time than the flow-angle method, especially if several 
streamlines are required as for a nozzle with a square test section. 

For the Mach number 10 nozzle small local irregularities in the 
streamlines were eliminated by integrating curves of tan 0 plotted 
against x where 0 was obtained from the stream- function method^ 

Foelsch (reference 7 ) has given analytic e3q)resslons for the 
coordinates of the transition streamlines in a three-dimensional nozzle. 
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Foelsch’s equations may be derived under somewhat more restrictive 
assumptions from equations (6) and (8) given in this paper. The 
points in the flow field eire Identified in the following sketch: 






R 


r 


cr 




\ (7+1)/2(7-1) 



(9) 


from eqviation (6). Then, by use of the relation 


^cr' sin (0 + p) 



cos 9 + 



1 sin 9 


and equation 


(8), equation ( 9 ) 


can be written as 
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where the integration is along the left Mach line PK. Making the 
assumption that M and 6 are constant along the left Mach lines and 
equal to their values at P and using the radial flow equations for 
rp/rcr and gives, for equation (lO), 


% ^P 


‘cr 


cr 


= J sin2 ^ 2( 


cos 


0p - co's ©cjlV^P - 1 sin 0p + cos 0pj 


1/2 


( 11 ) 


which is the same result as that given by Foelsch (equation (l7) of 
reference 7 ) • The x-coordinate of point K is then obtained by integra- 

d^B- 

V ^ 

tion of the relation d = along the left Mach line from 

rcr tan(0 + p) 

P to K on the assimption that 0 = 0p and M = as before. The 
result is 




— Sin Bp 


tan (0p + pp) 



cos 0p 


( 12 ) 


since 


eind 


^cr ^cr 


cos 0p 


Rr 


cr 


sin 0-n 

r 

•^cr 
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An approximation wMcli is more accurate for the Mach number 10 
nozzle than that given by Foe Is ch (reference 7 or equations (ll) and ( 12 ) 
herein) for the transition streamlines can be obtained from a considera- 
tion of the geometric properties of the flow in region IV (see sketch 
on p. 10 ) . Thus, if PK is assummed to be a straight line, (that is, 

M = Mp and 9=9^ along PK) , as was done in the derivation of 

equations (ll) and (12), the coordinates of point K are 




^P 

cos 0p 

Ask 


N 

^cr 

+ cos 

^cr 

\9p 

+ ^^pj 

rp 

sin 0p 

Ask 


\ 

^cr 

+ sin 1 

^cr 

[0p 

+ dpj 


(13) 


where ASj^ is the length of the Mach line PK. To obtain an expression 
for As, first consider its value for the final Mach cone DE. Thus, 


sin 



Re 

Ase 


ilh) 


If the flow were allowed to expand from point C as an undisturbed 
radial flow until it reached the spherical segment DD', the density and 
velocity would be the same at the surface of this spherical segment as 
it is along the final Mach cone DE. Therefore, conservation of mass 
requires that the surface area of the spherical segment DD' be equal to 
the plane cross-sectional eirea at E; that is. 


Be = 2rp sin ^ 


Then, from equation (lU), 


Ase = 2% rj) sin ^ 

The same reasoning may be applied to obtain an approximate relation 
between the area of the spherical segment formed by rotation of arc PP' 
about the radial line OP and the area Tt(KK')2, where the line KK' is 
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normal to the radial streamline extended through point P. Thus, 
assuming uniform flov along PK gives 


Bn — 0p 

ASk = 2Mprp sin — 


(15) 


This relation for Asg; can he obtained only hy assuming that the line 

OP is, in effect, the axis of a new nozzle with a final Mach cone of, 
half angle tip. This assumption is compatible with the supposition 
that the flow along PK is parallel and uniform. Substitution of equa- 
tion (15) into equations (13) results in . ' 


^ ^P 

^cr ^cr 

% = £L. 

^cr ^cr 
In equations (16) t 

flow in region III (fig. l) . As discussed previously this flow is 
determined by the specification of the final expansion angle and 

Mach number at point D. 

• The transition streamline through point C has been computed for the 
Mach number 10 nozzle by Foelsch's method (equations (ll) and (l2)) and 
the geometric method given by equations (16). The results of the two 
computations are con5)ared in table IV to the streamline coordinates 
obtained by the method of characteristics with the fine mesh. A com- 
parison between the streamline coordinates from the coarse and fine mesh 
is also included. The maximum error from Foelsch's method is about 
12 percent as compared to an error of h percent by the geometric method. 
Although both methods result in significant error for this nozzle, even 
in comparison with the values obtained by the characteristic method for 
the- coarse net, the possibility is that both methods may be more acciirate 
at lower Mach numbers where smaller values of 9q can be used. 


- 0p / \ 

cos 0p + Qlip sin cos (0p + 4pj 


^ (16) 


Bp — Pp 

sin Bp + 2Mp sin sin 


(®P + Pp) 


i values of rp, B-p, and Bq depend on the radial 


Cross-section contours.- The nozzle may be designed for a test 
section of arbitrary cross-sectional shape. Thus, the distance from the 
axis to a point on an arbitrary closed curve in a cixjss -sectional plane 
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of tlie final parallel and uniform flow region is the E ordinate of the 
starting point of a streamline which will he traced through the flow 
field to the minimum. Because of the symmetry of the flow field all 
streamlines hounding the nozzle can he designed in a single meridian 
plane. For the Mach number 10 nozzle a square cross section was 
specified in the uniform-flow region. Some of the resulting nozzle 
cross-section contours are shown in figure 4 . The cross section at the 
minimum is square since the flow at this station is assumed to he 
parallel and uniform. The maximum deviation of the other contours from 
a square is approximately 0.3. Nondlmensional coordinates for the 
Mach niimher 10 nozzle are presented in table V. The geometrical minimum 
is located at x^^ = 0.5043 r^^j. for reasons noted previously. 


CONCIUDIHG EEMAEKS 


A method is given for the design of three-dimensional nozzles 
utilizing axlsymmetric flow. The nozzle can he designed to produce 
parallel and uniform flow in a test chamber of arbitrary cross section. 
Design calculations are made for a Mach number 10 nozzle in which a 
square test section is specified to reduce the possibility of axlsym- 
metric Imperfections at the wall and to provide for the installation of 
schlieren windows. The specification of radial flow in a portion of the 
flow field in general assures a physically possible solution and results 
in an appreciable reduction of computing time. The flow in the Immediate 
vicinity of the axis of symmetry is computed by a simplified but suffi- 
ciently accurate method whereby the hodograph of any isentropic axlsym- 
metric flow is determined. The remainder of the supersonic flow field 
is constructed by the method of characteristics. The effect of mesh 
size on the accuracy of the calculation is obtained by comparing stream- 
lines throu^ two flow fields computed with steps of different size. 

This comparison Indicates for the Mach number 10 nozzle computation of 
this paper that the relative increase in accuracy is of the same order 
as the decrease in mesh size. Of the two methods considered for 
determining the final streamlines from the flow net, the stream -function 
method proved to be the most convenient and accurate. Foelsch's analytic 
expressions for the transition streamlines give coordinates which deviate 
from the flow net streamlines for the Mach munber 10 nozzle by as much 
as 12 percent. Similar analytic expressions derived in this paper from 
the geometric properties of the flow result in transition streamline 
coordinates with a maximum error of about 4 percent. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Via., April 1, 1952 
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appendix 

THREE-DIMENSIOKAL RADIAL FLOW 


For the steady flow of a gas in the absence of viscosity and heat 
transfer, there is obtained (reference 8, p. 553) 


• V = 1 f 
2 


(Al) 


if the flow along each streamline is assumed adiabatic. 


If the fluid moves radially and the velocity and speed of sound 
are functions of r only (where r is the distance from a fixed 
point), equation (Al) becomes 


(‘ 


Vr^yVr 



or, since Vj- is the only velocily component 



(A2) 


From the adiabatic energy equation the relation between Mach number and 
velocity is 


dV 1 dM 

V m2 M 

Then, combining eqviation (A2) and (A3) gives 


(A3) 


^ = i M^ - 1 
^ ' 2i + 2L^m2 m 


(A4) 



l6 
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Integratiiig eguatlon (Ak) 'between the limits r = 
M = 1 to M gives 


r _ 


= ^ 


L7 




^m2 


7+1 

2 ( 7 - 1 ) 


1 1/2 


■cr 


M 


to 


r and 


(A5) 


V- 


The characteristic egnations are now 
following sketch) : 


Left 


easily determined ( see the 
Mach line 


For the 



dx 


= tan(e + n) 


(A6) 


and 


dr 

r d0 


cot ^ 


Combining equations (Ak) and (A7) gives 


(A7) 


d0 


1 

2 


>/m^ - 


1 + m2 


dM 

M 


(A8) 


since cot n = i/m^ - 1. Integration of equation (a 8) from 9 = 0 to 0j 
and M = 1 to M gives 
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-(A9) 


Similarly, for the right Mach line. 


^ = tan(0 
dx 


m) 


(AlO) 


and 



i J7 . 

2 ^ 7-1 


tan 


-1 


7 + 


2(m2 - . 1 ) 


.-1I/2 


+ 1 


cos " — 


(All) 


Equations (A5) and (A9) have been computed and tabulated for 0j = 0° 
to 53.375° in intervals of 0.125° in table, I. Equation (A9) is used 
since the conventional definition of the total expansion angle is such 
that a positive increase in 'M results in a positive increase in ©j 
regardless of whether the Mach line is of the right or. left family. 

The Mach lines BC and CD (fig. l) may now be constructed as 
follows: The final Mach number is taken as = IO.068 which gives 

©I^ = 51*25° from table I. The flbw variables M and r/r^i* for 

any point on line CD are tabulated with the argument ©j = 

where the local flow angle © is the parameter. Thus, if © = i<-°, 

©X = i)-7.25°5 then from table I, M = 7-733272 and r/rcr = I2.77156. 

From this value of ^/^cr given value of © the coordinates 

of the point can be obtained since 


X _ r 
^cr ^cr 

B ^ r 
^cr ^cr 


cos © = 12.7^^ 


sin © = 0.8909 . 


For this particular nozzle, the flow angle at point C was taken as l6° 
which gives ©j^ = 35*25° and ©j^ = 19*25° since ©g =,0°. 
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TABIS I,- THREE 



2,557816 

543698 


25.20610 

25.10125 

22.99695 

22.89316 


2i.ga^ 

2 i.6ap5 

21.78^8 

21.68831 


21.50217 

21.20677 

21,11186 

21.01733 

20.92530 


2.860885 

2.87319I 

2l0979^ 

nm 


20.09552 

20.00557 


2.961040 

IMU 


19. 

19.47^ 
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TABIE II.- COMPARISON OF FLOW PARAMEirERS HEAR AXIS IN MACH 
NUMBER 10 NOZZLE AS COMPUTED BY EQUATIONS (3) AND (4) 


, Point , 
(fig. 2) 

Method 

x/^cr 

' cr 

M 

(deg) 

7,6 

Eq. (3) 
Eq. (4) 

Difference 

0.55477 

.55436 

.00041 

0.13186 

.12730 

.00456 

1.08526 

1.09874 

-.01348 

0.5889 

.59375 

-.00485 

6,5 

Eq. (3) 
Eq. (4) 

Difference 

.67197 

.67194 

.00003 

.08208 

.08153 

.00055 

1.23156 

1.23732 

-.00576 

1.0034 

1.0000 

.0034, 

5,1^ 

0 

Eq. (3) 
Eq. (4) 
Difference 

.82688 

.82686 

.00002 

.09904 

.09849 

.00055 

1.44344 

1.45212 

-.00868 

2.0000 

2.0000 

0 

S3 

Eq. (3) 
Eq. (4) 
Difference 

1.0352 

1.0352 

0 

.07293 

.07277 

.00016 

1.71835 

1.72354 

-.00519 

2.0000 

2.0000 

0 

3,2 

Eq. (3) 
Eq. (4) 

Difference 

1.2464 

1.2464 

0 

.06243 

.06234 

.00009 

2.00013 

2.00436 

-.00423 

2.0000 

2.0000 

0 

2,1 

Eq. (3) 
Eq. (4) 
Difference 

1.4737 

1.4737 

0 

.05731 

.05724 

.00007 

2.30450 
■ 2.30868 
-.004l8 

2.0000 

2.0000 

0 
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TABIE III.- COMPAZCSON OF ACCURACY OF THE FIOW 
IfETS ITT MACH MJMBER 10 NOZZLE FOR 0 = k^° 
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TABLE IV. - COMPARISON OF STREAMLIRE COORDINATES IN MACH 
NUMBER 10 NOZZLE FOR ^ = k'f 


^/^CT 


Vi-cr 

Characteristics method 

Geometric method 
(eq. 16) 

Foelsch method 


Coarse net 


Error 

(percent) 


Error 

(percent) 


Fine 

net 

Value 

Error 

(percent) 

Value 

Value 

3.7202 

1.067 

1.067 

0 

1.067 

0 

1.067 

0 

7.5 

2.016 

2.007 

-.45 

1.975 

-2.03 

1.919 

-4.81 

14.5 

3.280 

3.264 

-.49 

3.151 

-3.93 

2.956 

-9.88 

26,0 

4.624 

4.601 

-.50 

4.431 

-4.17 

4.067 

-12.05 

42.0 

5.706 

5.674 

-.56. 

5.523 

-3.21 

5.069 

-11.16 

57.0 

6.242 

6.197 

-.72 

6.121 

-1.94 

5.704 

-8.62 

73.5 

6.498 

6.469 

-.45 

6.457 

-.63 

6.203 

-4.54 

89.1219 

6.548 

6.548 

0 

6.548 

0 



6.548 

0 
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Flgvire 1 .- Scbsmatlc representation of the nozzle flov fields sliowinfe 
the Mach number variation with distance along the center line for 
the Mach number 10 nozzle . 
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Figure 2.- Mach net in region II for the Mach number 10 nozzle vith a 

streamline inc luded . 


mCA TR 2T11 








Figure 3-- Mach net in region IV for the Mach number 10 nozzle with a 

streamline included. 
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